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1 The Initial Results 

This note has its source in [4], in a different setting. Consider the system of 
recurrences, with ao = bo = 1, 

a n = On-i + kb n -i, (1) 

b n = On-l + K-l- ( 2 ) 

Later on we will generalize and show that the ratio tt is related to conver- 
gents to the square root of rationals. 

We postpone the proof by induction of the following formulas ("summation 
formulas" ) 
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Let a = k + 1 + 2\/fe, /3 = k + 1 - 2\/fc. Note that a + (3 = 2{k + 1), a/3 = 
(A; - l) 2 , a - (3 = 4\fk. Using Identities 1.87 and 1.95 in [2] and Pascal's 
Identity (see [3]) we have the following closed form representations 

a 2n = - + Vk , (7) 

a n+l +/? n+l 
«2n+l - ^ ' W 

n+l _ on+1 

^ = — *tf- (10) 

Let w(r, s) denote the recurrence 

w n = 2(k + l)w n -i - (k - l) 2 w n - 2 , w = r, wi = s. 

Furthermore let d n = w(l, k + 1), u n = w(0, 2k), v n = w(0, 2). The closed 
forms are 

_ a n + f3 n 
d n - ^ ' 

(a - P)(a n - (3 n ) 

U n - - , 

a n -I3 n 

V n = 2 — . 

a — p 

kv n = U n . (11) 

d2n = d n + U n , (12) 

d2n+l = d n +l, (13) 

b2n = d n + V n , (14) 

hn+l = V n+1 . (15) 

It follows that a>2n = w(l, 3k + 1), 02^+1 = w(k + 1, k 2 + 6k + 1), &2n = 
w(l, k + 3), b 2n +i = w{2, 4k + 4). Finally 

a n = 2(k + l)a n _ 2 - (k - l) 2 a n _ 4 



Note that 
Then 
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do = 1, ai = k + 1, d2 = 3fc + 1, 03 = k 2 + 6k + 1, 

6„ = 2(fc + l)6 n _ 2 - (fc - l) 2 6„- 4 

6o = 1, 6i = 2, 62 = + 3, 63 = 4fc + 4. 

These fourth order recurrences can be transformed in second order recur- 
rences in the following way. First of all note that 

l-2(k + l)x 2 + (k— 1) V = (l - 2x - {k - l)x 2 ) (l + 2x - (k - l)x 2 ) . 
The generating function of a n is 

1 + (fc + l)x + (k - l)x 2 - (fc - l) 2 x 3 
1 - 2{k + l)x 2 + (k - l) 2 x 4 ' 

This can be simplified to 

l + (fc-l)x 
l-2x-(k-l)x 2 ' 

so that we obtain 

a n = 2a n -i + (k - l)a n _ 2 , a = 1, ai = fc + 1. 

Analogously, the generating function of b n is 

1 + 2x - (fc - l)x 2 
\-2{k + l)x 2 + (fc-l) 2 x 4 ' 

and this becomes 
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1 - 2x - (fc - l)x 2 ' 

so that we obtain 

b n = 2b n -i + (fc - l)6 n _ 2 , 6 = 1, 61 = 2. 

Writing e = 1 + \/fc, 77 = 1 — \/fc (so that e = a?, rj = j3*) the closed forms 
are 

_ e n+1 + r] n+1 
a n - g ' 

,71+1 _ ^71+1 

T C '/ 

™ " 2^fc ' 



Now it is easy to show that lim n >00 = \fk. Indeed 



Or, 



e n+l + v n+l 



b n 2 e n+1 - 7? ri+1 

,71+1 I „n+l 

- V * 1 _ (2^+1- 

Since (2) n+1 converges to zero we have that ^ converges to y/k. 

Now we are ready to prove by induction the summation formulas. Assume 
that Equation 3, Equation 4, Equation 5 and Equation 6 hold for some 
integer n. This is tantamount to say that Equation 12, Equation 13, Equa- 
tion 14 and Equation 15 hold for some integer n. For n = they are satisfied 
since 

a = 1 = I fc u = d + u , 




a\ = k + 1 

6 = 1 = I : )/c u = do + v , 



Now 

a 2(n+l) = «2n+2 

= G2n+1 + kbin+\ 
= + &V n+ i 

= + U n+ i, 

where the third line is due to the induction hypothesis and we used Equa- 
tion 11. Then Equation 3 is satisfied through Equation 12. 

&2(ra+l) = ^2n+2 

= «2ra+l + &2n+l 

= d n+ i+v n+ i. 
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Again the third line is due to the induction hypothesis. Then Equation 5 is 
satisfied through Equation 14. 



a 2(n+l)+l = a 2n +3 

= (l2n+2 + &&2n+2 

= d n+ i + u n+ i + kd n+ i + kv n+ i 

= (1 + k)d n+ i + 2u n+ i 

a n+1 + {a n+1 - f5 n+l {a- 13) 

- (1 + k) - + - 

(a + (3)(a n+1 + (3 n+1 ) (a n+1 - /3 TO+1 (a - f3) 

2 + 4 

2a n+2 + 2(3 n+2 
I 

= d n+ 2- 

Here the third line is due to the fact that we already proved the formulas 
for a2n and &2n- Then Equation 4 is satisfied through Equation 13. For the 
last relationship consider first of all that 

(a n - P n )(a + (3) 

U n + V n = ~ 



a — (3 



Then 



&2(n+l)+l — &2n+3 



— d2n+2 + &2n+2 

= 2d n+ i + U n+ i + V n+1 

a n+i , pn+ i , (a n+1 -/3 n+1 )(a + )9) 

a — /3 

(a - /3)(a n+1 + /3 n+1 ) + (a n+1 - (3 n+1 ){a + (3) 



= 2 



a — j3 

a n+2 _ pn+2 



a-0 
= v n+2 . 

Here again the third line is due to the formulas already proven. Then Equa- 
tion 6 is satisfied through Equation 15. 
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2 Other Initial Conditions 



If we start out with ao = 0, bo = 1 and we denote the resulting sequences 
by a n and b n we get 



i=0 



i=0 



i=0 



i=0 



Then 



a n = 2a n _i + (A; - l)a n _ 2 , ao = 0, ai = fc, 
6 n = 26 n _i + (A; — l)6„- 2 , &o = 1, h = 1. 
It follows easily (n > 0) 

a n = kb n —i, 

bn = 0"n—li 

so that 2* — ► 

On 

3 Generalization 

The work done allows to analyze easily the following situation. Let us con- 
sider the sequences 

u n = u n -i + kv n -i, 

V n = hUn-i + V n -i, 

with k and h positive integers, where we start with no = 1, vq = 0. Then 
we get the following summation formulas and closed forms (where a\ = 
1 + hk + 2Vhk, p 1 = l + hk- 2^/hk) 

n ^nV 1 



u2n = J2[ 2i )(hky = -(a n 1+ ff), 



i=0 
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i=0 



-»= A g( 2i 2 :i)c*» i =^-^ 

»2„+i = ftf; ( 2 2 "| jj (hfc)' = - ft) + \ (a? + /ST). 

From this, using the same approach as before, we obtain 

u n = 2n„_i + (hk - l)u n - 2 , u = 1, u\ = 1, 

w n = 2w n _i + (/ifc - l)v n -2, v = 0, f i = /i. 
Using the closed forms given before we see that 

V 2n V ^ ' w 2n+l V h ' 

so that we can conclude 



4 Reduction 

Returning to the initial case, if k is an odd number the fraction ^ can 
be reduced and new sequences can be defined. More precisely, let k = 
2m + 1, m = 0, 1, 2, .... Then 

a = 2m + 1 + 1 + 2\fk 

= 2(771 + 1 + ^) 

= 2 7 , 

with 7 = m + 1 + Analogously, with 5 = m + 1 — Vk, we have 

/3 = 25. 

Note that 7 + 5 = 2(m + 1), 7< 5 = m 2 , 7 - 5 = 2\/fe = 2\/2m + 1. Now 
define 

/1\" 7™ + 5™ ^7" -5" 
c ^=[^j «2« = g + Vk 2 ' ' ' 
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C2n+1 = I g J a 2n+l = ^ , (17) 

„ 2 „ = (_) fe „ = LJ_ + 2_, (18) 
d 2n+ i = y Wi = 2 ^ • (19) 

Then of course 

d n 

Let u(r, s) denote the recurrence 

u n = 2(m + l)n n _i - m 2 u n - 2 , u Q = r, m = s. 

Then c 2n = 3m + 2), c 2 „+i = u(m + 1, to 2 + 4m + 2), d 2n = + 
2), d 2n+ i = u(l, 2(m + 1)). And finally 

c n = 2(m + l)c n _ 2 - m 2 c n _ 4 , 

d n = 2(m + l)d n _ 2 - m 2 d n -4, 

where the initial conditions are determined by the previous recurrences. 
The generating function of c n is 

2^3 



that of d n is 



1 + (1 + m)x + mx — m x 
1 — 2(1 + m)x 2 + m 2 x 4 



1 + x — mx 2 



1 — 2(1 + m)x 2 + m 2 x 4 

5 Some Identities 

We will work with the sequences 

a n = 2a n _i + {k - l)a n _ 2 , a = 1, oi = fc + 1, 

and 

6„ = 26 n _i + (fc - l)6 n _ 2 , b = l,h = 2, 
with closed forms 

£ n+l _|_ 
On = ^ i 
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,n+l _ ^n+l 

where e = 1 + y/k, r] = 1 — Using induction and Equations 1 and 2 we 
get 

a n = (k- l)6 n _i + 6„, 

q re+ i + (fc - l)a n -i _ , 

2A; " ° n ' 

fc6 n = a n + (k - l)a„_i. 

Using the closed forms (and the fact that e?? = 1 — k) we get 

c£-*# = (i-*) B+1 . 

a 2n = 2a n _ia n - (1 - k) n , 
(k — l)b m -ib n + b m b n+ i = b m+n+ i, 
(k — l)a m _ia n + a m a n+ i = kb m+n+ i. 
Writing n = m — 1 the last two become 

(k - + = b 2m , 

(k - l)a^_! + a 2 m = kb 2m . 
Using these identities we obtain 



2&a n /a 2 - (1 - k) n+l 2 
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°* = irfV^nr^ — rn- n (20) 

(1 - + 2fcfr 2 - 26 n ^kbl + (1 - fc)"+i 

&2n = ^ • (21) 

Now noting that a 2 n+i = 02.2™, writing 2 n instead of n we get recurrences 
for a 2 n and 62™ 



2fca2" /aL — (1 — k) 2 " +1 2aL ,„ ,,,» 



(1 - fc) 2 " +1 + 2fc6|„ - 2b 2 ^kb 2 2n + (1 - fc) 2n+1 
6 2 „ +1 = — 
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6 Newton's Iteration 

The Newton's iteration algorithm (see [8]) to approximate the square root 
of integers is given by 

1 / k 



starting with xq = 1. We assume k > 1. The limit of this recursion is the 
fixed point of the mapping 



which is \fk. 

For n > let us write 

_ On 

K 

We extend the sequences a n and b n setting cto = &o = 1. Then we have 

a n = a 2 n _ 1 + kb 2 n _ 1 , (22) 

b n = 2a n _i6 n _i. (23) 
We can generalize considering the recursion 

1 / k 



2 V hx n 

where h, k are positive integers with h ^ k. This produces an approximation 
to In this case we have 

a n = ha n _i + kb n _i, 

b n = 2ha n -\b n -\. 
We are going to prove using induction that, for n > 2, 

a n = 2a 2 n _ 1 -(k-l) 2n -\ (24) 

Using Equation 22 and Equation 23 we have a\ = k + 1, &i = 2, 02 = 
(fc+l) 2 +4fc = l+A; 2 +6fc. If we use Equation 24 we have a 2 = 2a\-(k-l) 2 = 
1 + k 2 + 6fc, so that Equation 24 is true for n = 2. Now assume that it holds 
for some n. Coupled with Equation 22 we get 

2a 2 _! - ttv^i = a 2 _! + fc6 2 _ l5 
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that is 

a n-l = kb l-l + w n-l, (25) 

where we wrote 

»n-l = (K - 1) 

Now 

= On + 

= ( 2 «n-i - ^n-i) 2 + 4&a 2 _ 1 & 2 _ 1 

= 4a^_! - 4a 2 _ 1 w n _i + w 2_ i + 4fc a 2_ 1 ^_ 1 . 

On the other hand 

2a 2 -w„ = 2{a 2 n _ 1 + kb 2 ^) 2 - w n 

= 2a 4 n _ 1 + 2k 2 b 4 n _ 1 + 4ka 2 l _ 1 b 2 n _ 1 -w n . 

Now using Equation 25 we have 

k 2 b n _ x = (a 2 ^ -u>„-i) 2 

= a^_ x + u> 2 _! - 2a 2 n _ l w n _ 1 . 

Then 

2a 2 - w n = 4a^_ x - 4a 2 _ 1 w„_i + 2w 2 l _ 1 + 4/ca 2 _ 1 6 2 _ 1 - w n 

since w ra = itf 2 _!. This concludes the proof. 
Using Equation 25 we get 



Cln-l = V^n-l + Wn-!, 

so that, through Equation 23, we get a recurrence for b n 

b n = 2b n ^ X \jkb 2 n _ 1 +w n - 1 . (26) 
Now we are going to prove, again by induction, the following closed forms 



a 2 " + 8 2 " 

a n = ±E , (27) 



or — 8 

6 " = ^v^ (28) 
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where a = 1 + Vk, ft = 1 — Vk. Note that a + (3 = 2, a/3 = 1 — k. Hence 
a n and b n are doubly exponential sequences (see [1]). 

For n = the closed form for ao gives = 1; for bo gives ^= = 1 so 
Equation 27 and Equation 28 are satisfied. Now assume that Equation 27 
and Equation 28 hold for some n. Then 

/ V"+/3 2 " \ 2 / g* n -p* n \ 2 

v 2 ; + v J 

a 2 " +1 + /3 2 " +1 + 2(a/?) 2 " a 2 " +1 + /3 2 " +1 - 2(a/3) 2 " 

"I - A; 



4 4fc 

cr + /r 



2 

This concludes the proof for a n . 
For 6 n we have 



b n +i — 2a n b n 



2Vk 

,2 n+1 o2 n+1 



a 1 _ - (3 1 
This concludes the proof. 

Incidentally we have proved, for n > 1, that kb n + w ra is a perfect square. 
With k = 2, a n is sequence A001601 and b n is sequence A051009 in [6]. 

Using Identities 1.87 and 1.95 in [2], where n is replaced by 2 n , we obtain 
the following summation formulas 



r=0 



on— 1 1 / 

2 n 



6 n = e ( 2r+ iK' n>0 - ( 3 °) 



From Equation 23 we see that 

n-l 

b n = 2 n J] Or, 

r=0 

which implies that 6 n is divisible by 2 n . 
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7 Related Sequences 

Let us consider the sequence c n , with co = l, c± = r, r > 1, such that for 
n > 2 

c n = 2c„_ 1 — 1. 

From [5] we get the closed form 

7 2 " + S 2n 
c n +\ — , 



where 7 = r + \/r 2 — 1, 5 = r — y/r 2 — 1. Define the sequence <i„ by 

7 2 -<5 2 

«n+l 



2Vr 2 - 1 ' 

We set do = 1 and we obtain d\ = 1. Note that d 2 = 2r. Then, for n > 1, 

on — 1 „on — 1 on — 1 r on-l 

7 +5 7 -<5 

«n+l 



2 2^ 2T7 T 

2c n c? n . 



Then 

n-l 

dn = 2™- 1 n a. 
i=l 

From this we can evaluate 



lim n( i + - 



n >oo 

i=l 



Indeed we have 



n , 1 



i=l 



1 + Cl 1 + C2 1 + C n 

C\ C 2 C n 

1 + ci 2cl 2c 2 _! 

Cl c 2 c n 
(l + ci)2 ra - 1 cic 2 ---c re -i 

(r + 
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Using the closed forms of c n and d n it is easy to see that 

lim — ^ = y/r 2 — 1. 

Hence 

lim nfi + 1 )=^ = ^- 

The case r = 3 is considered in [7]. 
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